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Abstract. We give a constructive proof that coordinate transformations ex¬ 
ist which raise the regularity of the gravitational metric tensor from (7°’^ to 
in a neighborhood of points of shock wave collision in General Relativ¬ 
ity. The proof applies to collisions between shock waves coming from different 
characteristic families, in spherically symmetric spacetimes. Our result here 
implies that spacetime is locally inertial and corrects an error in our earlier 
RSPA-publication, which led us to the false conclusion that such coordinate 
transformations, which smooth the metric to C^’^, cannot exist. Thus, our re¬ 
sult implies that regularity singularities, (a type of mild singularity introduced 
in our RSPA-paper), do not exist at points of interacting shock waves from 
different families in spherically symmetric spacetimes. Our result generalizes 
Israel’s celebrated 1966 paper to the case of such shock wave interactions but 
our proof strategy differs fundamentally from that used by Israel and is an 
extension of the strategy outlined in our original RSPA-publication. Whether 
regularity singularities exist in more complicated shock wave solutions of the 
Einstein Euler equations remains open. 


1. Introduction 

The guiding principle in Albert Einstein’s pursuit of General Relativity (GR) 
was the principle that spacetime should be locally inertial [T], (we say also locally 
Minkowski or locally flat). That is, an observer in freefall through a gravitational 
field should observe all of the physics of Special Relativity, except for the second 
order acceleration effects due to spacetime curvature (gravity). But the assump¬ 
tion that spacetime is locally inertial is equivalent to assuming the gravitational 
metric tensor g is smooth enough so that one can pursue the construction of 
Riemann Normal Goordinates at a point p\ he., coordinates in which g is ex¬ 
actly the Minkowski metric at p, such that all first order derivatives of g vanish 
at p, and all second order derivatives of g are bounded in a neighborhood of p. 
However, the Einstein equations are a system of partial differential equations 
(PDE’s) for the metric tensor g coupled to the sources and the Einstein equa¬ 
tions by themselves determine the smoothness of the gravitational metric tensor 
by the evolution they impose. Thus the condition on spacetime that it be locally 
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inertial at every point cannot be assumed at the start, but must be determined 
by regularity theorems for the Einstein equations. 

The presence of shock waves makes this issue all the more interesting for the 
Einstein equations with a perfect fluid source. In this case the Einstein equations 
G = kT imply the GR compressible Euler equations DivT = 0 through the 
Bianchi identities, [2], and the compressible Euler equations create shock waves 
whenever the flow is sufficiently compressive, [3]. At a shock wave, the fluid 
density, pressure, velocity, and hence T are discontinuous, so that the Einstein 
equations imply the curvature G must also become discontinuous at shocks. But 
discontinuous curvature by itself is not inconsistent with the assumption that 
spacetime be locally inertial. For example, if the gravitational metric tensor 
were (7^’^, (differentiable with Lipschitz continuous hrst derivatives), then second 
derivatives of the metric are at worst discontinuous, and the metric has enough 
smoothness for there to exist coordinate transformations which transform g to 
the Minkowski metric at p, with zero derivatives at p, and bounded second 
derivatives as well, |1]. Furthermore, Israel’s theorem asserts that a spacetime 
metric of regularity (7°’^, (i.e., Lipschitz continuous), across a smooth single 
shock surface, is lifted to (7^’^ by the (7^’^ coordinate map to Gaussian normal 
coordinates, and this is smooth enough to ensure the existence of locally inertial 
coordinate frames at each point [S]. 

In [B], Groah and Temple set out a framework in which to address these issues 
rigorously by providing the hrst general existence theory for spherically symmet¬ 
ric shock wave solutions of the Einstein-Euler equations allowing for arbitrary 
numbers of interacting shock waves of arbitrary strength. In coordinates where 
their analysis is feasible. Standard Schwarzschild Goordinates (SSG)j3the gravi¬ 
tational metric is only (7°’^ at shock waves, and it has remained an open problem 
as to whether the general weak solutions constructed by Groah and Temple could 
be smoothed to (7^’^ by coordinate transformation, as was proven by Israel for 
single shock surfaces, [5]. 

In this paper we partially resolve the open problem of Groah and Temple by 
proving there do exist coordinate transformations that lift the regularity 
of the gravitational metric tensor from (7°’^ to G^’^ at a point of shock wave 
interaction between shocks from different characteristic families in spherically 
symmetric spacetimes. In [7] the authors introduced the idea of a regularity sin¬ 
gularity, a point in spacetime where the metric tensor is (7°’^ but not G^ regular 
in any coordinate system. Our result here is the hrst step in extending Israel’s 
theorem to interacting shock waves, by proving that spacetime is indeed locally 
inertial and that no regularity singularity exists at points of such shock wave col¬ 
lision. This negates our false conclusion in [7] that regularity singularities exist 
at points of shock interaction, and the error in [7] is explained and corrected in 
this paper, (c.f. Lemma ITT]) . The question as to whether regularity singularities 
can be created in more complicated solutions of the Einstein Euler equations, by 
more complicated shock wave interactions, remains an open problem. 

spherically symmetric metric can generically be transformed to SSC, c.f. [5] 
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The existence of regularity singularities would be surprising for General Rel¬ 
ativity, where it is commonly assumed that the gravitational metric tensor is at 
least The metric regularity is the threshold regularity required for the 
existence of locally inertial coordinate frames, the existence of which is essential 
for properties of shock waves in Minkowski space-time to be recovered in the limit 
of weak gravitational helds. The metric regularity is a starting assumption 
in the singularity theorems of Hawking and Penrose, [8]. At a regularity singu¬ 
larity the metric would be Lipschitz continuous but not in any coordinate 
system, so discontinuities in the metric derivatives would be present in every 
coordinate system, and this would open the door for possible new gravitational 
effects. The authors will address the implications of regularity singularities in a 
forthcoming paper. 

To state our main result precisely, let denote a spherically symmetric 
spacetime metric in SSC, where the metric takes the form 

ds^ = g^ydx^dx^' = —A{t, r)dt^ + B{t, r)dP -|- r'^dVP'. (1.1) 

At the start either t or r can be taken to be timelike, and dVt^ = -f-sin^('d)(i(/9^ 

is the line element on the unit 2-sphere, c.f. [6j. In Section [3l we make precise 
the dehnition of a point of regular shock wave interaction in SSC between shocks 
from different families. Essentially, this is a point in {t, r)-space where two shock 
waves enter and leave a point p, such that the metric is Lipschitz continuous 
across the shocks and smooth away from them, the Rankine-Hugoniont (RH) 
jump conditions hold across each shock curve and are continuous up to the point 
of interaction p, derivatives of all quantities are continuous up to the shock 
boundaries, and the SSC Einstein equations hold weakly in a neighborhood of 
p and strongly away from the shocks, |9]. The main result of the paper is the 
following theorem: 

Theorem 1.1. Suppose that p is a point of regular shock wave interaction be¬ 
tween shocks from different families, in the sense that conditions (i) - (iv) of 
Deftnition \3.1\ hold, for an SSC metric g^y. Then the following are equivalent: 
(i) There exists a C^'^ coordinate transformation in the {t,r)-plane, 

with Jacobian Jjf, defined in a neighborhood Af of p, such that the metric 
components ga/s = JfJgd^iv cire C^'^ functions of the coordinates x". 

(a) The Rankine Hugoniot conditions, fl3.4p - 03.51) . hold across each shock 
curve in the sense of (v) of Deftnition \3.1[ 

Furthermore, the above equivalence also holds for the full atlas of C^'^ coordinate 
transformations, not restricted to the {t,r)-plane. 

Our proof of Theorem 11.11 provides an explicit method for constructing the 
Jacobians of {t, r)-coordinate transformations that smooth the components of 
the gravitational metric from (7°’^ to in a neighborhood of p. In order to 
prove the equivalence in Theorem 11.11 we characterize all such Jacobians that 
lift the metric regularity. Our method of proof differs substantially from the one 
used by Israel, the latter being based on studying the Einstein tensor in Gaussian 
normal coordinates and concluding that the metric is C^'^ in these coordinates. 
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The main ideas for the proof of Theorem 11.11 were already introduced in [7], but 
an error in the last section led us to the wrong conclusion that metric-smoothing 
is not possible. In fact, Sections [2] - O of this paper mostly coincide with the 
corresponding sections in [7] . In Section [TJ we correct the error in [7] and outline 
the proof of Theorem 11.11 

Our assumptions in Theorem 11.11 apply to shock wave interactions in which two 
timelike shock waves enter and leave the point of interaction. In the case of a per¬ 
fect fluid source, this type of interaction is realized between two incoming shock 
waves from different characteristic families, c.f. |9]. Although points of shock 
wave interaction are straightforward to construct for the relativistic compressible 
Euler equations in flat spacetime, and a general existence theory for shock wave 
interactions in GR is given in [B] , we know of no complete mathematically rigor¬ 
ous construction of a point of shock wave interaction in GR sufficient to derive 
detailed structure at such points. However, our assumptions regarding regular 
shock wave interaction in SSG given in Dehnition 13.II below, are straightforward, 
are consistent with [6], and confirmed by the numerical simulations in |10j . 

In Section [3] we set out the framework of shock waves in GR, and define 
what we call a point of regular shock wave interaction in SSG. In Section H] we 
introduce a canonical form for functions (7°’^ across a hypersurface and write the 
RH conditions as a relation between first order metric derivatives. In Section O 
we derive necessary and sufficient conditions that Jacobians of coordinate 
transformations lift the regularity of the metric tensor from to at points 
on a single shock surface. This is the so-called smoothing condition. In Section 
[6] we give a new constructive proof of Israel’s theorem for spherically symmetric 
spacetimes, combining the results from Sections H] and |5l For this, we hrst derive 
a canonical form of the Jacobian which satisfies the smoothing condition and 
show that the freedom to add an arbitrary G^-function to our canonical form 
suffices for the Jacobians to be integrable to coordinate transformations. Sections 
2-6 here agree with Sections 2 - 6 in [7], in principal. 

In Section [7] we extend the above strategy to the case of shock interaction 
and outline the proof of Theorem 11.11 The content of this section is new and 
complete details can be found in Section 7 - 10 of m- For the proof, we first 
derive the canonical form for the Jacobian satisfying the smoothing condition 
across each of the shock curves. This is the source of the error in our original 
RSPA paper [7], since we incorrectly omitted terms which encode the presence 
of two shocks. As in Section [6l we now need to prove that one can integrate 
the Jacobian to coordinates. This is achieved by showing that the freedom 
to add a G^ function to the canonical form of the Jacobian suffices to solve 
the integrability condition, ^ which is done as follows: choosing two 

of the free functions arbitrarily, (say <I>^ and the integrability condition 
turns into a linear first order system of PDF’s of the form Ut + cUr = F{U) 
for the remaining two free functions U = ($q,<I>q) as unknowns. The source 
term F{U) is non-local, depends on the restriction of U to the shock curves 
and derivatives of U along the shocks, and is discontinuous at the shock waves. 
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We first prove existence of a (7°’^ solution U, which by itself is not a sufficient 
smoothness for the resulting Jacobian to meet the smoothing condition. We 
then use a bootstrapping argument to show that U is indeed regular if and 
only if the RH jump conditions hold. (Interestingly, there is an apparent loss 
of smoothness across the characteristic which passes through the point of shock 
wave interaction, even though this characteristic curve lies within the region of 
smoothness of the SSC metric. The resulting metric is but seems to be no 
smoother. This is a kind of mild new irregularity by itself and further indicates 
the subtlety of the problem.) The above construction gives us Jacobians which 
smooth the metric before and after the interaction takes place, but to obtain 
the Jacobian in a spacetime-neighborhood of the point of interaction, we need to 
match these two Jacobians across the surface t = 0, the time of shock collision. 
Thus, in the last step we prove that one can choose the free functions at f = 0 
appropriately for the metric in the resulting coordinates to match across the 
f = 0 interface such that it maintains its regularity. Interestingly, again the 
RH conditions come in at this hnal step to ensure that this matching can be 
done consistently. In this construction of the Jacobian several conditions, which 
appear over-determined at the start, are consistent as a consequence of the RH 
jump conditions, giving us conhdence that all terms have now been accounted 
for. 


2. Preliminaries 

Let g denote a Lorentzian metric of signature (—1,1,1,1) on a four dimensional 
spacetime manifold M. We call M a C^-manifold if it is endowed with a C^- 
atlas, a collection of four dimensional local diffeomorphisms from M to such 
that any composition of two local diffeomorphisms x and y of the form x o y~^ is 
regular. The mapping x o y~^ is referred to as a coordinate transformation. 
In this paper we consider C'^’^-manifoIds. 

Our index notation for tensors use Greek letters p, z/,... G {t, r, 9, (p} for SSC 
coordinates (in which the spacetime metric g is (7°’^) and Greek letters a, /?,... G 
{0,1, 2, 3} for transformed coordinates, gap = We use the Einstein 

summation convention whereby repeated up-down indices are summed over all 
values for the given indices. Tensors transform by contraction with the Jacobian 
= 1^, the inverse Jacobian is denoted by J", and indices are raised and 
lowered with the metric and its inverse g^^, which transform as bilinear forms, 
g^iv = J^Judap- We use the fact that a matrix of functions is the Jacobian of 
a regular local coordinate transformation if and only if the curls vanish, i.e., 

= and (2.1) 

where f^a = ■§^ denotes partial differentiation with respect to the coordinate x" 
and Det{JI^) denotes the determinant of the Jacobian, c.f. [2]. 

In this paper, we do not restrict to, but are motivated by the Einstein Euler 
equations, 




( 2 . 2 ) 
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which couples the metric tensor to the undifferentiated perfect fluid sources 

= (p + p)u^^u'^ + (2.3) 

through the second order Einstein curvature tensor where 

div T = 0 (2.4) 


follows from div G = 0. Here k = SttQ, where Q is Newton’s gravitational 
constant, p is the energy density, Ui the 4-velocity, and p the pressure, c.f. 00 
Equation fl2.4p reduces to the relativistic compressible Euler equations when 
is the Minkowski metric, and the Euler equations close when an equation of 
state (e.g. p = p{p)) is imposed. Shock waves form from smooth solutions of 
the relativistic compressible Euler equations when the initial data is sufficiently 
compressive, [9]. 

Across a smooth shock surface S, the RH jump conditions hold. 


[rnn. = 0 , 


(2.5) 


where [/] = fi — fn denotes the jump in / from right to left across S, and rij, is 
the surface normal. The RH condition fl2.5p is equivalent to the weak formulation 
of fl2.4p across S, c.f. [S]. 

In this paper we restrict to time dependent spherically symmetric metrics in 
Standard Schwarzschild Coordinates where the metric takes the form fll.ip . The 
Einstein equations for a metric in SSC are given by (c.f. 0 ) 


with 


B^ + B 


- A- 


B - 1 
r 

Bt 
B- 1 


Bit — Ar 


$ 


-kAB'^tT^^ 
nAB^rT^^ 
-2KABr^T^^ , 


BAtBt B^ Ar , ABr , A^ ^ ArBr 
2AB ~ ^ 74 ^ 2B 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 


Note that the first three Einstein equations in SSC imply that the metric cannot 
be any smoother than Lipschitz continuous if the source T is discontinuous (for 
example, G T°°), and in this paper we make the assumption throughout 
that A and B are Lipschitz continuous, i.e., functions of t and r. 


^The Riemann curvature tensor introduced in differs from the one used by us and in [8] 
by a factor of —1 which, in [2], is compensated for by setting k = —SttQ. MAPLE uses the sign 
convention in [2] for the Riemann tensor, which is important to keep in mind when computing 
the Einstein tensor for (12.6p - (|2.9I) with MAPLE. 
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3. A POINT OF REGULAR SHOCK WAVE INTERACTION IN SSC BETWEEN 
SHOCKS FROM DIFFERENT FAMILIES 

In this paper we restrict attention to radial shock waves, by which we mean 
hypersnrfaces S locally parameterized by 

S(t,d,(p) = (3.1) 

across which A and B are and T in fl2.3l] satishes 02.51) . Then, for each t, 
S is a 2-sphere with radins x{t) and center r = 0. Treating cj) and 9 as constant, 
we introdnce 7 , the restriction of a shock snrface S to the (f, r)-plane, 

'lit) = (3.2) 

with normal 1 -form 

n„ = {x,-l). (3.3) 

For radial shock snrfaces 03.11) in SSC, the RH jnmp conditions 02.51) take the 
simplihed form 

[T°°] X = [T°^] , (3.4) 

[Ti°] X = [T^^] . (3.5) 

To generalize the above framework to collisions between shocks from different 
families, we think of the incoming and ontgoing branches of the two shock waves 
as fonr distinct timelike shock snrfaces, parameterized in SSC by 

d, (j)) = (t, xf{t), 9, 0), (3.6) 

with i = 1,2, and where is dehned for f < 0 and S)*" for f > 0 . Assnme Jlf 
intersect at f = 0 , that is, 

3^^(0) = To = x^(0), 

for some tq > 0. Restricted to the (f,r)-plane, are described by the shock 
cnrves 

jfit) = {t,xf{t)), (3.7) 

with normal 1 -forms 

inf)u = (3.8) 

We assnme the 'yf are with all derivatives extending to f = 0. Denoting with 
[•])*^ the jnmp across the shock cnrve the RH conditions now read. 


A“]fV = 

AT. 

(3.9) 

[ryfif = 

AT- 

(3.10) 


For the proof of Theorem 11.11 it snffices to restrict attention to the lower 
(f < 0) or npper (t > 0) part of a shock wave interaction that occnrs at f = 0. 
That is, it snffices to consider the lower or npper half plane in separately, 

= {(f, r) : f < 0} or R^ = {(t, r) : f > 0} , (3-11) 
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respectively. (We denote with the closure of M|.) Whenever it is clear that 
we restrict consideration to or we drop the superscript ± of the quantities 
introduced in fl3.6p - fl3.10p . 

We now dehne “a point of regular shock wave interaction in SSC between 
shocks from different families” as a point p where two shock waves collide, re¬ 
sulting in two outgoing shock waves, such that the metric is smooth away from 
the shock curves and Lipschitz continuous across each shock, allowing for a dis¬ 
continuous and the RH condition to hold. In the special case that 
describes a perfect fluid, this type of collision corresponds to an interaction of 
shock waves in different characteristic families, c.f. [S]. 

Definition 3.1. Let tq > 0, and assume to he an SSC metric in (7°’^ (A/"), 
where M <Z ^ is a neighborhood of the point p = (0,ro) of intersection of the 
timelike shock curves i = 1,2, introduced in fl3.7p . Let J\ denote the open set 
consisting of all points in M not in the image of any . Then, we say that p is 
a “point of regular shock wave interaction in SSC between shocks from different 
families” if: 

(i) The pair {g, T) is a strong solution of the SSC Einstein equations fl2.6p - fl2.9l) 
in Af, with T^^'-' e C^{J\f) and g^u G C^{Af). 

(a) The limits of T^^^ and of the metric derivatives g^v^a exist on both sides of 
each shock curve 'yf, including the point p. 

(in) The jumps in the metric derivatives [gfj,^,a]fit) are C^ functions for all 
fG(—e, 0] or for all t G [0,e). 

(iv) The (upper/lower)-limits 

exist. The (upper/lower)-limits exist for all derivatives of [gi^uAf- 

(v) The stress tensor T is bounded on Af and satisfies the RH conditions 

irntiaX = 0 

at each point on 'yf{t), t G (—e, 0) or t E (0,e), and the limits of these 
jumps exist up to p as t ^ D. 

The framework introduced above mostly agrees with our original setting in 
[7] . However, in contrast to the above definition, in [7] we imposed the structure 
only on or R/ separately, because there we looked for a contradiction, while 
here we look for a construction. 

4. Functions across a hypersurface 

In this section we give a precise dehnition of functions that are C^’^ across a 
hypersurface and use this to derive a canonical form for such functions. 

Definition 4.1. Let E be a smooth (timelike) hypersurface in some open set 
Af C We call a function f “Lipschitz continuous across T”, (or C^’^ across 
fl), if f G C^'^{Af), f is smooth (f G C‘^{Af\fl) suffices) inAf\fl, and limits of 
derivatives of f exist and are smooth functions on each side of S separately. We 
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call a metric Lipschitz continuous across S in coordinates if all metric 
components are across S. 

The main point of the above definition is that we assume smoothness of / 
away and tangential to the hypersurface S. Note that the continuity of / across 
S implies the continuity of all derivatives of / tangent to S, i.e., 

[f,aV = 0, (4.1) 

for all tangent to S. Moreover, Dehnition 14.11 allows for the normal derivative 
of / to be discontinuous, that is, 

[f,aK ^ 0, (4.2) 

where is normal to S with respect to some (Lorentz-) metric dehned on 

M. 

We can now clarify the connections between the Einstein equations and the 
RH jump conditions fl3.4p . fl3.5p for SSC metrics that are only across a 
hypersurface. To this end, consider a spherically symmetric spacetime metric 
fll.ip given in SSC, assume that the first three Einstein equations fl2.6p - fl2.8p 
hold, and assume that the stress tensor T is discontinuous across a smooth radial 
shock surface described in the (f,r)-plane by 7 (f) as in fl3.ip - fl3.3p . Condition 
fl4.ip across 7 applied to each metric component g^^, in SSC fll.ip then reads 


[Bt] = -x[Bri (4.3) 

[At] = -x[Ar\. (4.4) 

On the other hand, the hrst three Einstein equations in SSC fl2.6p - fl2.8p imply 

[Br] = (4.5) 

[Bt] = (4.6) 

= KAB'^r[T^^]. (4.7) 


Now, using the jumps in Einstein equations fl4.5p - fl4.7p . we find that fl4.3p is 
equivalent to the first RH jump condition fl3.4p . (c.f. Lemma 9, page 286, of [1]), 
while the second condition fl4.4p is independent of equations fl4.5p - fl4.7p . because 
At does not appear in the first order SSC equations fl 2 . 6 p - fl 2 . 8 p . The result, then, 
is that in addition to the assumption that the metric be across the shock 
surface in SSC, the RH conditions fl3.4p and fl3.5p together with the Einstein 
equations fl4.5p - fl4.7p . yield only one additional condition over and above fl4.3p 
and fl4.4p . namely, 

[A] = -x[Bt] . (4.8) 

The RH jump conditions together with the Einstein equations will enter our 
method in Sections |MZ1 only through the three equations fl4.8p . fl4.3p and fl4.4p . 

The following lemma provides a canonical form for any function / that is 
Lipschitz continuous across a single shock curve 7 in the (f, r)-plane, under the 
assumption that the vector normal to 7 , is obtained by raising the index in 
fl3.3p with respect to a Lorentzian metric g that is (7°’^ across 7 . (Note that by 
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Definition 14.11 varies in directions tangent to 7 . Here, we suppress the 
angular coordinates.) 

Lemma 4 . 2 . Suppose f is across a smooth curve j{t) = {t,x{t)) in the 
sense of Definition \4.1l t G (—e, e), in an open subset M o/M^. Then there 
exists a function <l> G C^{Af) such that 

r) = r), (4.9) 

where 

(oJMUL ec\-t,e). ( 4 . 10 ) 

n^Ufj 

and n^{t) = {x{t),—l) is a 1-form normal to the tangent vector v'^{t) = j^{t). In 
particular, it suffices that indices are raised and lowered by a Lorentzian metric 
which is across 7 . 

In words, the canonical form fl4.9p separates off the kink of / across 7 (that 
is, the (7°’^ element of /) from its more regular behavior away from 7 : The 
kink is incorporated into |T(f) — r\, (p gives the smoothly varying strength of the 
jump, and <h encodes the remaining behavior of /. 

In Section [7] we need a canonical form analogous to fl4.9p for two shock curves, 
but such that it allows for the Jacobian to be in the weaker regularity class (7°’^ 
away from the shock curves. To this end, suppose two timelike shock surfaces 
described in the (t,r)-plane by, ji{t), such that fl3.6p - fl3.10p applies. To cover the 
generic case of shock wave interaction, we assume each ji(t) is smooth (at least 
(7^) away from t = 0 and all derivatives extend continuously to f = 0. It suffices 
to restrict to upper shock wave interactions in 

Lemma 4 . 3 . Let ji{t) = {t,Xi{t)) be two smooth curves defined on I = (0,e), 
for some e > 0, such that fl3.6p - fl3.8p hold. Let Af be an open neighborhood 
of P = (0,ro) in and suppose f is in fl M^), but such that f is 

tangential to each 7 * with (14.ip holding. Then there exists a (7°’^ function 4) 
defined on A/" fl such that 

[•hi], = 0 = [$,]„ (4.11) 

for i = 1,2, and 

/(^i - r| + ^{t, r), (4.12) 

i=l,2 

for all {t, r) in M where 

( 4 . 13 ) 

and {ni)^{t) = (xj(f),—1 ) is the 1-form normal to vf{t) = ji(t), for i = 1,2, 
and indices are raised by a Lorentzian metric (7°4 across each %. 
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5. A Necessary and Sufficient Condition for Smoothing Metrics 


In this section we derive a necessary and sufficient pointwise condition on 
the Jacobians of a coordinate transformation that it lift the regularity of a 
metric tensor to in a neighborhood of a point on a single shock surface S. 
This is the starting point for our methods in Sections [6] and [71 Proofs and further 
results can be found in Section 5 of 
We begin with the transformation law 


9a0 



(5.1) 


for the metric components at a point on a hypersurface S for a general 
coordinate transformation x“, where, as customary, the indices indicate 

the coordinate system. denotes the Jacobian of the transformation, that 
is, = 1^. Assume now, that the metric components are only Lipschitz 
continuous with respect to across S. Then differentiating fIS.ip in the direction 
w = we obtain 




(5.2) 


where [/] denotes the jump in the quantity / across the shock surface S. Thus, 
since both g and are in general Lipschitz continuous across S, the jumps 
appear only on the derivatives. Equation fl5.2p gives a necessary and sufficient 
condition for the metric g to be in x" coordinates. Namely, taking w = 
in SSC, fl5.2p implies that [ga 0 ,'y] = 0 for every a,/J, 7 = 0 ,..., 3 if and only if 


= 0 . (5.3) 

Note that if the coordinate transformation is (7^, so that is (7^, then the 
jumps in J vanish, and fl5.2p reduces to 


which is tensorial because the non-tensorial terms cancel out in the jump [ga 0 ,'y]- 
It is precisely the lack of covariance in fl5.2p for (7^’^ transformations that provides 
the necessary degrees of freedom in the jumps [J^^] to lift the smoothness of a 
Lipschitz metric one order at a single shock surface. 

We now exploit linearity in fl5.3p to solve for the [7^ o-] associated with a given 
(7^’^ coordinate transformation. To this end, suppose we are given a single radial 
shock surface S in SSC locally parameterized by 




(5.4) 


For such a hypersurface in Standard Schwarzschild Coordinates (SSC), the angu¬ 
lar variables play a passive role, and the essential issue regarding smoothing the 
metric components by (7^’^ coordinate transformation, lies in the atlas of (f, r)- 
coordinate transformations. Thus we restrict to the atlas of {t, r)-coordinate 
transformations for a general (7°’^ metric in SSC, c.f. fll.lj) . The following lemma 
gives the unique solution of fl5.3p for (f, r)-transformations of (7°’^ metrics 
g in SSC. 
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(5.5) 


Lemma 5.1. Let 

= —A{t, r)dt^ + B{t, r)dr‘^ + r‘^dVt^ , 

he a given metric expressed in SSC, let S denote a single radial shock surface 
fl5.4p across which g is only Lipschitz continuous. Then the unique solution [J^„] 
of fl5.3p which satisfies the integrability condition in 5'5'C'Jl c.f. fl 2 .ip . 

W.M = [JW 

is given by: 

[4,t] = 

[Jit] = 

W = 

[Jit] = 


A 

A 

[^r] 

B 

[^r] 

B 


Jn + 


Jl + 


Jr. + 


A 

A 

K 

B 


■ r 




r 

•^0 


0 ,.] 


Jt [^jr 

Jl+ 


[Jl 

[Jl, 

[Jl 

[Jl 


A 

[^ 

A 

B 

B 


Jn + 


& 

A 

[St 


jr 

'^0 


r -L ^ ■' r 


Jn + 


A 

B 

[S: 


r 


Jl + '-^Jl ) -(5.6) 


{We use the notation /i, p, cr e {t,r} and a,(3 E {0,1}, so that t,r are used to 
denote indices whenever they appear on the Jacobian J.) 


To avoid confusion in Section [H] and [71 we introduce the notation 

XL = 1A,J (5-7) 

to denote the right hand sides in fIS.Op . 

Condition fl5.3p is a necessary and sufficient condition for [gap,-I = 0 at a 
point on a smooth single shock surface. Because Lemma 15.11 tells us that we 
can uniquely solve fl5.3p for the Jacobian derivatives, it follows that a necessary 
and sufficient condition for [gap,'f\ = 0 is also that the jumps in the Jacobian 
derivatives be exactly the functions of the jumps in the original SSC metric 
components recorded in fl5.6l) . Thus, Lemma [5.11 implies the following lemma: 


Lemma 5.2. Let p he a point on a single smooth shock curve 7 , and let g^y 
he a metric tensor in SSC, which is (7°’^ across 7 in the sense of Definition 
\4.1\ Suppose is the Jacobian of a coordinate transformation defined on a 
neighborhood Af of p. Then the metric in the new coordinates gap is in C^'^{N') 
if and only if J^f satisfies fl5.6p . 


6. Metric Smoothing on Single Shock Surfaces and a 
Constructive Proof of Israel’s Theorem 


In this section we outline an alternative constructive proof of Israel’s Theorem 
for spherically symmetric spacetimes, (see Section 6 in m for complete details). 
For the proof, in light of Lemma [5^ we need to construct Jacobians of coordinate 
transformations, defined in a neighborhood of a point on a single shock surface, 

^We use here that = J^ is equivalent to (EH), but with derivatives taken in 

SSC, c.f. [TT] for more details. 
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that satisfy fl5.6p . In other words, we need to introduce a set of functions, J^, 
that satisfies fl5.6p and the integrability condition fl2.ip in some neighborhood of 
the shock. The main theorem of this section is the following: 

Theorem 6.1. (Israel’s Theorem) Suppose g^i, is an SSC metric that is 
across a radial shock surface S in the sense of Definition \4- 1\ such that it solves 
the Einstein equations fl2.6p - fl2.9p strongly away from E for a which is 
continuous away from E. Let p be a point on E. Then the following is equivalent: 

(i) There exists a coordinate transformation of the {t,r)-plane, defined in 
some neighborhood M of p, such that the transformed metric components 
are functions of the new coordinates. 

(a) The RH conditions, fl3.4p - fld.Sp . hold on EH J\f' for some J\f' D M. 

Furthermore, the above equivalence also holds for the full atlas of coordinate 
transformations, not restricted to the {t,r)-plane. 


The main step is to construct Jacobians acting on the (f, r)-plane that satisfy 
the smoothing condition fl5.6p on the shock curve, the condition that guarantees 
[Ooig^'y] = 0. The following lemma gives an explicit formula for functions 
satisfying fl5.6p . The main point is that, in the case of single shock curves, both 
the RH jump conditions and the Einstein equations are necessary and sufficient 
for such functions to exist. 


Lemma 6.2. Let J\f be a neighborhood of a point p, forp lying on a single shock 
curve 7 across which the SSC metric g^y is Lipschitz continuous in the sense 
of Definition \4.1\ and let g^y be defined on Af. Then, there exists functions 


G which satisfy the smoothing condition fl5.6p on 7 fl A/" z/ and only 

if the RH conditions fl4.8p hold on jHAf. Furthermore, any such function is 
of the “canonical form” 


Jaifc) = IT(f)-r|+ <!)(((t,r) (6.1) 


with 

<(() = (6-2) 

where is defined in fl5.7p . /i G {t,r}, a G {0,1}, and G C'°’^(A/') satisfy 

[a,<i>((] = o = [a,«i>;(]. (6.3) 

Explicitly, the Jacobian coefficients are given by 


Toit) 

T\{t) 

Toit) 

T[it) 


[Ar](j){t) + [Bt]u{t) 
AA o 7(t) 

[Ar]iz{t) + [Bt\C{t) 
AA o 7(t) 

[Bt](j){t) + [Br]uj{t) 
AB o 7(t) 

[BtHt) + [Br]at) 

AB o 7(t) 


(6.4) 
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where 

0 = $* O 7 , a; = 'I'o O 7 , // = $* o 7 , ( = $!; o 7 . (6.5) 

Furthermore, the above equivalence also holds for the full atlas of coordinate 

transformations, not restricted to the {t,r)-plane. 


Proof. Suppose there exists a set of (7°’^ functions satisfying fl5.6p . It is shown 
in m that these functions satisfy 


[A,l = -iK.r] 

for all n G {t,r} and a G {0,1}. Combining fl6.6p for the special case jj, 
a = 0 with the right hand side in fl5.6p leads to 


1 

2 



[^r] jr 


X 

2 



A 



( 6 . 6 ) 
t and 


Using now the jump relations for the metric tensor, fl4.3p - fl4.4p . hnally gives 
[Ar] = —x[Bt], that is, the non-trivial RH condition fl4.8p . 

For proving the opposite direction, as a consequence of Lemma 15.11 it suffices 
to show that all t- and r-derivatives of the functions Jjf, dehned in (16.11) . satisfy 
05.61) for all p G {f, r} and a G {0,1}. Observing that 06.5p implies the identities 

0 = J* O 7, 1/ = o 7, w = o 7, ( = J[ o 7 , 

and using the matching of the functions 06.3p . as well as the RH condi¬ 

tions in the form 04.3p . 04.41) and 04.81) . it follows immediately that the Jacobian 
ansatz 06.11) satishes 05.6p . This proves the existence of functions Jjf satisfying 
the smoothing condition 05.61) . Finally, applying Lemma [4.21 it follows that all 
functions satisfying 05. 6 p assume the canonical form 06.ip . 

In HH, Section 6, we extend the Lemma beyond coordinate transformations 
in the {t, r)-plane and the corresponding canonical form is described. □ 


To complete the proof of Israel’s Theorem, we need to show that there exist 
functions such that the J^, dehned in (16.ip . satisfy the integrability condition, 
02. ip . that is, p = Jp a- consider and as given functions 

and write (12.ip as a PDF in the unknown 


and it is straightforward to show that 02.ip is equivalent to the system of PDF’s 
dtU + cdrU - MU = {\X\M - H{X) (h - c) ) ^ - |X| ^ , (6.7) 


where X{t, r) 


x{t) — r and the coefficients are given by 


c = 


JJ 


= ^ and 




P 

jr 


P 

P 

"' 1 ,! 


The goal now is to prove we can solve 06.7p for U G C^{J\f) D C^{Af \ 7 ). 

Fquation 06.71) is a system of non-local PDF’s, since the right hand side of 
06.71) contains the Jacobian coefficients Jq and ip^ which depend on U o 7 itself. 
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and standard existence theory cannot be applied right away. Nevertheless, pre¬ 
scribing initial data on the shock curve, the right hand side of fl6.7p becomes a 
given source term and 06.71) turns into a strictly hyperbolic linear system of hrst 
order PDE’s. Imposing the condition 

c 7 ^ XV, (6.8) 

which ensures that the shock curve is non-characteristic, the standard existence 
theory in [12] proves existence of a solution U of 06.71) . by integration of the initial 
data and the source term along characteristic lines, c.f. m, Section 6 for more 
details. The existence theory yields a solution U which lies in C°’^(A/')nC'^(A/'\ 7 ) 
and is smooth away from 7 , but it does not give us the necessary regularity 
across the shock, 06.3p . due to the presence of the (discontinuous) Heaviside 
functions H{X) in 06.7p . The hnal step to complete the proof of Israel’s Theorem 
is now to show that the coefficients of H{X) in 06.7p vanish on the shock curve 
precisely when the RH jump conditions hold, as stated in the next lemma, which 
then yields the desired regularity across 7 . 

Lemma 6.3. Assume the assumptions of Theorem \6.1\ and denote with f and h 
the coefficient functions of the Heaviside function H{X) in the first and second 
component of 06.7p . respectively. Then, 

/o 7 = 0 = /io 7 (6.9) 

if and only if the RH conditions, 03.4p - 03.Sp . hold on 7 . 

Proof. To derive an explicit expression for the coefficients to H{X) in 06.71) . note 
that the matrix AA contains Heaviside functions as well. Then, collecting all 
terms containing H{X) and using X o 7 = 0 and 06.Sp . we hnd 

/07 = LplxV-Lp\xf) + ip\u-Lpl(, 

ho'y = ip^x V — ip'lx (f) + (fluj — iPq(. (6.10) 

Now, replace ifh and (p\ by their dehnition, 06.4p . then a straightforward compu¬ 
tation shows that / o 7 = 0 is equivalent to 

([H^] x[Bt\) {(fC - = 0. (6.11) 

Now, using 

{( pC - vuj ) = det { J ^ o ^) ^ 0 , (6.12) 

we conclude that / o 7 = 0 if and only if [A^] -|- x[Bf\ = 0, which is equivalent 
to the RH condition 03.5|) . c.f. 04.81) . Similarly, replacing and in O6.10p by 
06.41) . a straightforward computation shows that /i o 7 = 0 is equivalent to 

{[Bt] + x[Br]) (0C - = 0. (6.13) 

Now, using again 06.121) . the equivalence of h o 7 = 0 and 04.4|) follows, and thus 

to the RH condition 03.4p . This completes the proof. □ 

We can now complete the proof of Israel’s Theorem. As shown above, there 
exist functions such that the dehned in 06.11) solve the integrability con¬ 
dition, 02.ip . Moreover, by Lemma 16.31 these have the required regularity. 
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that is, they satisfy fl6.3p at the shock curve and away from the shocks. By 
Lemma 06.21) . the satisfy the smoothing condition 05.6p if and only if the RH 
jump conditions hold. Taken all together, we constructed a Jacobian which is 
integrable to coordinate functions and which maps the Lipschitz continuous SSC 
metric to a metric that is C^’^regular in the resulting coordinates if and only if 
the RH jump conditions hold. (See |TT], Section 6 , for the complete proof.) This 
proves Theorem [ 6 l 

7. Metric Smoothing around Points of Shock Wave Interaction 

In this section, we outline the proof of Theorem 11.11 the details of which can 
be found in nn, Section 7-10. In principal, we follow the ideas from the 
previous section: We hrst extend our Jacobian ansatz fl 6 .ip to the case of two 
interacting shock waves and then show that this set of functions can be integrated 
to coordinates. However, in contrast to the single shock case addressed in the 
previous section, we have to pursue the construction of on the upper and 
lower half-plane, separately, and then show that the resulting functions can 
be “glued” together in a way appropriate to smooth the metric. 

We now begin constructing the Jacobian. In contrast to the single shock case, 
the restriction ° Ji, i = 1 , 2 , does not only depend on the free functions 
but also on the Jacobian coefficients from the other shock curve, j ^ z, that is, 

°7* = (7-1) 

The error in [7] was to falsely omit the term |3:i(-)—X 2 (-) | in fl7.1l) . Since the 
smoothing conditions fl5.6p depend on o itself, and since {<^2)0 depends on 
and {(pi)a on ( 7 ’ 2 )a, we have to prove that the Jacobian coefficients, 
are well defined, in the sense that they can be consistently defined in terms of 
the metric and the free functions alone. The following proposition is the key 
step in extending Israel’s Theorem to shock interactions. It gives the canonical 
form of Jacobians that meet the smoothing conditions across each shock curve 
in either or and act on the (f,r)-plane only. Without loss of generality 
we formulate the proposition for 

Proposition 7.1. Let p be a point of regular shock wave interaction in SSC 
between shocks from different families, in the sense of Definition 1,9. Jl with (i) - 
(iv) being met, with corresponding SSC metric, defined on ^17®+. Then 
the following are eguivalent: 

(i) There exists functions Jf G (7°’^ (^Af fl M+j, for p G {t, r} and a G {0,1}, 

which satisfy the smoothing condition fl5.6p on 7 * fl A/" fl M+, for i = 1 , 2 . 

(a) The RH condition fll.Sp holds on each shock curve 74 ( 7 ^ 7 ®+, for i = 1,2, 
as in Definition \3.1l (v). 

Furthermore, any such set of functions Jf is of the “canonical form” 

Jaifr) = \xi{t)-r\ + ^>f{t,r), 

i=l,2 


(7.2) 
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where G yj\f fl j have matching derivatives across each shock curve 
liif), for t > 0 , that is, 

[dr^>f], = 0 = [dtK]i G Va G {0,1}, (7.3) 

and where (<y 2 i)a defined implicitly through 

(»>(): = (7.4) 

with {Ji)ar denoting the values in (15.Oh with respect to 7 *. Explicitly, the 
values for are given by 


{^1% = (7-5) 
{^iyi = (7-6) 


i^^Yo = 


iv^Y = 


+ [Br]j 


1 - BijBji 


i5l(|B.li<l>‘ili + [Brl.>I>;li) + 45-([ST'I>‘U + lBrh<S>’l\,)B, 


1 _ S B- 


(7.7) 


(7.8) 


with j ^ i in i\7.7\\ and (I7.8p . and where we define Ai = A o -j^, Bi = B o 
and 


(7.9) 




|(ri(-) - 3:2(-)| 
4R; 


[B,|i - [B,|. I . 


(7.10) 


Furthermore, the above eguivalence also holds for the full atlas of coordinate 
transformations, not restricted to the {t,r)-plane. 


The proof of Proposition 17.11 is recorded in [TT], Section 7. To reiterate, the 
error in [7], was to falsely omit the term ((Pj)^ |ti(-) — T 2 (-)| in 117.ip . which we 
correct here. The effect on the Jacobian coefficients fl7.7p - fl7.8p is precisely the 
appearance of the non-zero function Bij, and fl7.7p - fl7.8p reduce to the (incorrect) 
formulas in [7] upon setting Bij = 0. 

Following the argument in Section [ 6 l the next step in the construction is to 
first write the integrability condition as a PDF in the unknown U = ^(<hQ,<FQ), 
considering d)} and as given smooth functions which enter the coefficients of 
the PDF, and then to prove existence of a suitable regular solution U. As in 
Section [ 6 l we write the integrability conditions as 

dtU + cdrU = F{U), 


c 


Jl 

JV 


M 


L( Jit Jir\ 

•7} V Jit •^l,r ) 


where 


(7.11) 
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and, setting Xi{t, r) = Xi(t) — r, 

F = HiX,) (i. - c)) ( (^1 ) - IX.I ( ) } . 

Again, the difficulty proving the existence of solutions to fl7.1ip is that F{U) 
contains the non-local source terms ((pi)o and ((pi)o, which depend on f/ o yj 
and its derivatives. In contrast to the single shock case, localizing the F{U) by 
imposing initial data on the shock curves is problematic because of the lack of 
regularity at the point of interaction and the subsequent gluing problem. For this 
reason, we develop an iterative scheme in which, we replace F{U^) by F{U^~^), 
where is the known prior iterate at the k-th step. One of the problems 

proving convergence of this scheme is controlling the derivatives in F{U’^~^), 
which are of leading order in fl7.1ip . These terms can be controlled, since all 
these derivatives are multiplied by \xi{t) — r|, z = 1,2, which are small close to 
the point of interaction. In Section 8 of we prove that these iterates indeed 
converge to a Lipschitz solution of fl7.1ip . In generalization of Lemma 16.31 it 
is then shown that this solution has the crucial regularity across the shocks, 
fl7.3p . necessary for the construction of the Jacobian smoothing the metric tensor 
from (7°’^ to c.f. Proposition 17.11 

In Section 8 of m we prove that these iterates indeed converge to a solution 
of fl7.1ip with the required regularity. This is accomplished by hrst proving 
the existence of a Lipschitz continuous solution and then bootstrapping to the 
crucial regularity. This bootstrapping requires the RH jump conditions, as 
in the proof of Lemma 16.31 but a further difficulty is that the regularity along 
the charactersitic curve emanating from the point of shock interaction must also 
be addressed. See Section 8 in HU for details. The result is recorded in the 
following proposition, we formulate it in but the same result holds in 

Proposition 7.2. Assume regular initial data t/o(r) and assume and 
are given functions. Then, there exist a neighborhood M of p and there exist 
a regular function U = ^(<Fq,<Fq) which solves the integrability condition, 
fl7.1ip . in the region such that U{0,r) = Uo{r) for all (0,r) G M+. 

Proposition 17.21 dualizes the construction of Jacobians on which are in- 
tegrable to coordinate transformations and map the Lipschitz continuous SSC 
metric to a regular metric in the new coordinates. A similar construction 
gives Jacobians with the same properties on To complete the proof of The¬ 
orem 11.11 it remains only to prove that one can “glue” these Jacobians at the 
{t = 0)-interface and maintain the metric regularity, by a suitable choice of 
the free functions at t = 0. 

For this, we hrst introduce some notation. We denote all objects in fl7.2p - 
fl7.10p with an additional index “-I-” or ” to indicate whether they are dehned 
in or respectively. For instance, denotes the canonical Jacobian on 
A/± = AAnMl, {<pf)a its coefficients and its free functions. We denote with 
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{■} the jump across the (t = 0)-interface, that is, 

juKr) = limuft, r) — limuft, r), for r ^ rn, 
t/'O t\o 

{ii}(ro) = um- - UM+-, (7.12) 

where u is some function for which the above limits are well-dehned, e.g., the 
metric or the Jacobian derivatives, and where um- denotes the limit at p of 
u restricted to the region in between the two shock curves and um+ the 
respective limit between the upper two shock curves. 

We now derive the conditions for matching conditions which are necessary 
and sufficient for the metric in the new coordinates, gap = to be 

regular across the (t = 0)-interface. The condition that the Jacobian matches 
continuously across the {t = 0)-interface is 

= ./r(O.r) - JU(0,r) = 0, (7.13) 

for all r. We call these the -matching conditions, c.f. [TT] . 

For the matching of the Jacobian derivatives, we follow the reasoning in Sec¬ 
tion E] which leads to the smoothing condition, fIS.jp . but now we apply this 
reasoning to the {t = 0)-interface. That is, the condition that gap is continu¬ 
ously differentiable across the {t = 0)-interface is given by 

{gapA = 0. (7.14) 

Substituting gap = fl7.14p and using fl7.13p as well as the SSC- 

metric being regular away from the shocks, i.e., {fi'^y,cr}(r) = 0 for all r ^ tq, 
we conclude that the -matching conditions are given by 

({A,}A + {yU}A)9p. = 0, Vrjlro, (7.15) 

+ atr = r„. (7.16) 

To finish the proof of Theorem 11.11 we now show that fl7.13l) and fl7.15p - 
fl7.16p are met for the Jacobian in fl7.2l) . by appropriately matching the free 
functions and as well as their t-derivatives at t = 0. At the start, the 
conditions (17.1311 and (17.1511 - (I7.16P appear over-determined, essentially because 
the derivatives dtLh^ = are not free to assign, but determined 

by equation (17.1 ip at t = 0. Nevertheless, (17.1111 together with the RH jump 
conditions give the consistency of the matching conditions, (I7.13P and (17.1511 - 
(I7.16p . within the freedom available and imply exactly the matching of and 
This is all achieved in the final lemma, the proof of which is recorded in 
[TT] . Section 9. 

Lemma 7.3. Let be two Jacobians of the canonical form (17.211 . defined on 
A/± = A7nR^ respectively, with corresponding free functions Assume that 
the integrability condition ( 12 .ip holds and that Jj( 0 ,r) 7 ^ 0 , by an appropriate 
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choice of $i(0, r). If the match att = t), such that 

KKr) = 0, (7.17) 

{S,4S}(r) = -((vjr)‘i + (%)'i-W)‘i-W)‘i)|r-r„|, (7.18) 

{a,4>;}(r) = - ((v^r);+ (%);-(tip;-( 4 iL;)|r-r„| (7.19) 

hold for all ( 0 ,r) G J\f, then satisfies 07.131) and fl7.15p - 07.16p . 

Since we are free to choose the initial data of 07.lip as well as the remaining 
free functions, and $ 1 , such that 07.17P - 07.19P and J((0, r) 7 ^ 0 hold, Lemma 
17.31 implies that our canonical Jacobian can be matched such that the metric in 
the new coordinates is regular. This completes the proof of Theorem 11.11 

8. Conclusion 

Our result shows that no regularity singularities exist at points of shock wave 
interaction between shocks from different characteristic families, and this corrects 
the false conclusion in [7]. We prove that one can extend Israel’s result to shock 
wave solutions of the Einstein equations containing points of shock wave interac¬ 
tion between shocks from different characteristic families in SSC. We introduce a 
new method for constructing Jacobians of coordinate transformations that map 
the SSC metric to a regular metric in the new coordinates. Our method 
differs fundamentally from Israel’s proof based on Gaussian normal coordinates. 
Israel’s proof does not extend because Gaussian normal coordinates do not ex¬ 
ist in a neighborhood of a point of shock wave interaction. Whether regularity 
singularities can be created by shock interactions more complicated than the 
interaction of two spherical shock waves in SSG from different families, remains 
an open problem, even assuming spherical symmetry. But our method opens 
the door to address this regularity issue for more complicated solutions. Perfect 
fluid matter models are essential to the description of astrophysical phenomena. 
Thus the question as to the existence of regularity singularities is fundamental 
to General Relativity, both because their existence would change the mathemat¬ 
ical framework for GR perfect fluids, and because they would give rise to new 
detectable astrophysical effects. 
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